Density Perturbations in the Brans-Dicke 

Theory 

J. P. Baptista 1 , J. C. Fabris 2 and S. V. B. Gongalves 3 
Departamento de Ffsica 
Universidade Federal do Espfrito Santo 
Goiabeiras - CEP 29060 -900 
Vitoria - Espfrito Santo 
Brazil 

February 7, 2008 



1 E-mail: plinio@cce.ufes.br 
2 e-mail: fabris@cce.ufes. br 
3 E-mail: sergio@cce.ufes.br 



Abstract 



We give here the calculation of density perturbations in a gravitation theory 
with a scalar field non-minimally coupled to gravity, i.e., the Brans-Dicke 
Theory of gravitation. The purpose is to show the influence of this scalar 
field on the dynamic behaviour of density perturbations along the eras where 
the equation of state for the matter can be put under the form p = ap, where 
a is a constant. We analyse the asymptotic behaviour of this perturbations 
for the cases a = 0, a = —1, a — 1/3 and p = 0. In general, we obtain a 
decaying and growing modes. In the very important case of inflation, a — — 1, 
there is no density perturbation, as it is well known. In the vacuum phase the 
perturbations on the scalar field and the gravitational field present growing 
modes at the beginning of the expansion and decaying modes at the end of 
this phase. In the case a = it is possible, for some negative values of uj, to 
have an amplification of the perturbations with a superluminal expansion of 
the scale factor. We can also obtain strong growing modes for the density 
contrast for the case where there is a contraction phase which can have 
physical interest in some primordial era. 



1 Introduction 

The existence of a classical scalar field in Nature has been considered in many 
theories of gravitation that present alternatives to General Relativity. The 
prototype of scalar theories is the Brans-Dicke Theory 0, 0,0,0, whose 
Lagrangian is given by: 



It can be expected that the presence of the scalar field leads to different 
predictions with respect to those we obtain in General Relativity. In cos- 
mology, this can lead to far reaching consequences, since the standard sce- 
nario given by General Relativity presents, besides some spectacular success, 
important drawbacks such as the horizon, flatness and structure formation 
problems[|J. 

However, local physics limits the value of the parameter u: it must be 
greater than 500 to account the classical tests. This result has reduced the 
interest in the Brans-Dicke Theory. This situation has changed recently 
with the proposal of the extended inflation in the de Sitter phase, 

Brans-Dicke Theory predicts power-law inflation instead of exponential; in 
fact, following D. La and P. Steinhardt, at the beginning of inflation the BD 
solutions (for p = —p) approaches the Einstein-de Sitter solution. In the 
second stage of the inflation, both the scalar field and the scale factor grow 
by power law rather than exponential. This feature prevents the so-called 




(1) 
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"graceful exit" problem. However, in order to work, the parameter must be 
u) m 24, contradicting observation; nevertheless this constraint follows from 
local conditions. This drawback can be overcome through a generalization of 
the original Brans-Dicke Theory, allowing the parameter uj to be a function 
of the field <fi itself. 

This revival of Brans-Dicke Theory leads us to ask if the problem of struc- 
ture formation can be modified through the introduction of scalar field. This 
question has been treated in many different situations in the literature [||,[|], 
nj. Here, we propose to study the evolution of density fluctuation in the 



traditional Brans-Dicke Theory in the different phases of the Universe. Even 
if we consider wasa constant, this analyses can furnish many insights in how 
the scalar field modifies the main conclusions about gravitational instability 
obtained employing General Relativity. In particular we will see that in de 
Sitter phase, the scenario differs substantially from the traditional one. In 
the other phases, however, the differences are much less important. As it was 
mentioned above, the main interest in the presence of the scalar field is the 
possibility that this field accelerates the growth of density perturbations. 



We will work in the Lifschitz-Khalatnikov formalism [11]. The reason is 
that it furnishes essentially the same relevant results we can obtain employing 
another formalism, and also it allows for the choice of synchronous frames 
where the physical interpretation of the concerned quantities are easily done. 
We must, of course, be careful about the presence of unphysical modes [|12 



However they can be eliminated by performing an infinitesimal coordinate 
transformation. Taking care of the so called residual coordinate freedom, we 
can be sure to retain just the physical modes. 



We will alllow for negative values for the parameter u. In fact, a negative 
to is what is predicted by the effective models comming from Kaluza-Klein 
and Superstring theories [p~3j| . 

This paper is organized as follows. In section II we present the back- 
ground solutions of the unperturbed Universe in the following phases of its 
development: vacuum, inflation, radiation and incoherent matter. In sec- 
tion III we obtain the perturbed equations and their solutions in terms of 
Bessel functions for that phases. In section IV we calculate the asymptotic 
behaviour of the solutions for t — > and t — > oo. Finally, in section V 
we discuss the relations between the wavelength A and the particle-horizon 
distance H~ x . 

We will assume in this article the following notations: the greek indices 
run from zero to three; the latin indices run from one to three; the signature 
is (+,—,—,—); we use a Robert son- Walker metric with flat space section 
(k = 0); the scalar field is a time function; the energy-momentum is the 
perfect fluid; we use the synchronous gauge which fix the reference frame. 

2 Background Solutions 

We assume that background Universe is spatially flat, homogeneous and 
isotropic, i. e., it is described by the Robertson- Walker metric: 

ds 2 = dt 2 - a 2 dx 2 , (2) 
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where the a(t) is a scale factor of the Universe, and c = 1. The energy- 
momentum tensor of the background matter takes a perfect fluid form: 

= (p + pWU" -pg" v , (3) 

with an equation of state p = ap. In these expressions p is a density matter, 
p is a pressure and a is a constant. 

From the Lagrangian density ([I]), we obtain the field equations: 

+ t(</w-#^ d 0) , (4) 

□0 = — — — T . (5) 

Inserting the metric (Q), into equations (f4|) and (|5|) we obtain the equa- 
tions of motion: 

'd 8ir [2 + uo + 3a + 3auj\ uo - 2 1 • ■ . N 

- 3 a = J P { 3T^ J V + / ' (6) 

+ 3-0=-^-p(l-3a) . (7) 

The above equations must be supplemented by a conservation equation: 

T^ ]v =Q ; (8) 

when p = we have the energy conservation: 

p + 3-p(l + a) = , (9) 
a 

while, for p = j we get the condition: 

p,z = . (10) 
4 



Background solutions can be obtained assuming that the scale factor a(t) 
and the scalar field <fi(t) have a power-law form: 

a{t) oc f , <j)(t) oc t s . (11) 

By direct substitution of the above relation in the equations (|7|), (^) 
and (0), w e have: 

9 „ /2 + toil + 3a) + 3a\ 9/ . . 

- 3r(r - 1) = Svrpt 2 - 5 J +s\l +u)-s , (12) 

s 2 + s(3r _ 1) = _|!L_ pt 2- s(1 _3 a) ; (13) 

p + 3^p(l + a) = , (14) 

p,t = . (15) 

In what follows we shall consider the important special cases: vacuum, 
inflation, radiation and incoherent matter. 



2.1 Vacuum (p = 0) 

Here, the Universe has no ordinary matter and the energy-momentum tensor 
is null in all space-time. Our background solutions from equations fljjD , ( [13]) 
are 

S = l-3r , r= 1 ^ , (16) 

where u > — |. 
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2.2 Inflation (a = -1) 

In this scenario we have a drastic expansion of the Universe during the early 
period of the Big Bang. In this case, the background solutions are: 

s = 2 , r = u+ l - . (17) 

Strictly speaking, inflation occurs for (jJ~>\. 

2.3 Radiation (a = |) 

In this era the energy density of the Universe was tottaly dominated by 
relativistic particles. The solutions are: 

8 = , T = \ ■ ( 18 ) 

As it is easily seen the above solutions are the same for the General Relativity 
by the same scenario ||14|| . 



2.4 Incoherent Matter (a = 0) 

This phase is such that the Universe is dominated by nonrelativistic matter 
with negligible pressure. The solutions are: 

2 + 2uo 

s = 2 - 3r , r = . 19 

3 Perturbed Equations 

Our objective here is to calculate the perturbations of the equations (|), (0) 
and (S), and to obtain the exact solutions for the perturbed equations for 



6 



each scenario of the Universe. We are interested in the density matter per- 
turbations so that we will only use the temporal components of the Einstein's 
equations since that component is directly associated with the scalar modes. 
Adopting the synchronous gauge, we put h ^ = 0. 

In order to derive the perturbed equations, we write the metric tensor as: 



(20) 



leading to the Ricci tensor 



(21) 





(22) 



the energy-momentum tensor trace is 



T = T + 5T = T + 5p-35p ; 



(23) 



for the d'Alembertian of the scalar field we find 



50(f) = 5(f> + aah' 



<fi - —h kk <j) + 3-5<fi -5(f) . 



(24) 



cr a a 



Using now the following definitions 

i) h kk = a 2 h , 

ii) 5(f> = X(f> with A << 1 , 
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iii) 5p = Ap with A << 1 , 

we will get the set of perturbed equations: 

In A 



(3 + 2u>) 



p(l-3a) , (26) 



A-(l + a)[h-5U k , k 1=0 , (27) 



1 + «)(2-3a)-5[/ i + 5U j (l + a) = -a^aAj/ . (28) 



a 
a 

Equations fl25|), (p5[) are the differential equation for the perturbation 
h(t) of the gravitational field, X(t) is the scalar field perturbation and A(t), 
the perturbation for ordinary matter. Equations (|27|) and ( p8|) are the per- 
turbed equations derived from the conservation equation (§). With the above 
set of differential second order equations, we will analyse each phase of the 
Universe. 

In what follows, we will suppose that the perturbations have a plane wave 
behaviour: 

X(x, t) — A(i)exp(— iq.x) and h(x, t) = h(t)exp(— iq.x) , (29) 
where q is the wavenumber of the perturbations. 

3.1 Vacuum (p = 0) 

In the Universe without matter there is only perturbations of the metric and 
scalar field. The two coupled linear second-order equation for h(t) and X(t) 



are given by: 



^h + ~h = \ + 2-(l + u){l - 3r) 



1 • 



1 



h(l - 3r) = At + 2A(1 - 3r) + 3Ar + -5— [<?A . 



(30) 
(31) 



2 -v- , - ■ v > ' ,2 

Differentiating (|31|) and removing /i and /i, we find the following linear 
third-order differential equation for A(i): 

"2 



A +A 



t 



l + r) 



+ 



+A 



£(16 + 12w) - £(24 + 18w)^(2 + 2a;) + ^ 



9' 



(32) 



In order to solve the above equation, we perform the following substitu- 
tion: 



A 



K 
1 



k = t r 7 , x = t p , with p = 1 — r , 



7 = x m g , with m 



We get from equation ([32]), 



1 - 2r 
1 — r 



yY + yg' + 9(y 2 -i) = o 



xq 
1 — r 



(33) 



i.e., Bessel equation of order 1, where the prime denotes derivatives with 
respect to the conformal time. 

The solutions of (|33"D are the Bessel and Neumann functions of the first 
kind: 

g = C 1 J 1 (y)+C 2 N 1 (y) . (34) 
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The solutions for X(t) and h(t) are given by: 



h(t) 



+ 



[1 - 3r)t 2 
1 

(1 - 3r)t r 
+ 



t 1_r (Ci^i(y) + C 2 N 1 (y)^dt + C 3 



Qr + 



2q 2 



f 2r-4 



4 — Qr — -(1 + r 



CiJi{y) + C 2 N x (y) 



(35) 



+ 



+ 



C 1 j 1 (y)+C 2 N 1 (y) 



dt + C. 



4 j 



(36) 



1 - 3r)F 

where Ci,C 2 ,C 3 and C 4 are constants and Jk,N k means respectively the 
time derivative of Bessel and Neumann functions. 



3.2 Inflation (a = -1) 

The remarkable feature of this case is that the density contrast is null. This 
can be seen through equations (|2T|)-(|28|). So, we have two coupled linear 
second-order equation for h(t) and X(t): 

2r ■ -4 • lGn 

h+ / l _2A--(2r+l)A--— -rp\ = , (37) 
t t (1 + r) 

2 + 3r\ . 1 2x 1; 16tt 



x + 2X \—) + (^ qX -l h + WT7f x -°- (38) 

The resolution method is very similar to the vacuum case. So we find the 
following linear third-order differential equation for X(t): 
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By employing the transformations: 



K 



A = — , k — f*7 , x = t p with p = 1 — r , 

V 



7 = x m g with m = 



1 - 7r 



2(1 - r) ' ^ 1-r ' 
we obtain the Bessel equation of the order k: 



xq 



y 2 g" + yg' + (y 2 -k 2 )g = , 



where 



k = 



r + 3 



2(1 -r) • 
Its solution is 

g{y) = CMy) + C 2 J- k (y) , 

where Jk(y) is the Bessel function of the order k. 
The solutions for X(t) and h{t) are given by: 

1 



X(t) 



t 



h(t) 
(3r + 2) 



( Ci J fe (y) + C 2 J_ fe (y) )dt + C 3 
1 -5r 



l + 5r 



(d J fc (|/) + C 2 J_ k (y))- 



1 — 5r 



+ 



- \ ( / * ^ + C 2 J- k {y))dt + C 3 ) - 



1 

+ 7 



+ 



(6r + 2) 1 2 
t 2 ^t 2 ^ 



*— (0x^(2/) + C 2 J_ fc (y)) 
| t^ L (C 1 J k (y)+C 2 J-M)dt + C 3 



where Ci, C 2 , C 3 , C4 are constants. 
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In contrast with the vacuum case, r is greater than 1 for uj > ~. For 
r = 1 {uj = |), equation (|39"D becomes an Euler equation and admit the 
particular power-law solution: 



X = t n , 



n = -2 ± a/4 - q 



(44) 



3.3 Radiation (a = |) 

By considering the background solutions for the perturbed quantities, we 
have, from equations 



(45) 

(46) 
(47) 
(48) 



A + §A + ^A = , 

2 • 4 
A = -h--5U\ J , 



5U j + 2t5U j 



1 



A ,3 



Taking equation ( |4B] ) and performing the following transformations: 
x = t p , with P = 2 ' ^ = 2 ' m ^ w ith m = — - and ?/ = 2ga; 



we obtain the differential Bessel equation: 



y y 2 



9 = 



(49) 



The solution for \(t) is: 

A(t)=H(c 1 J i (y) + C7 a J_j(i/)) 
where Ci and C*2 are constants. 



(50) 
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Now, from equations p5|),p7|) and the equation (f|8|),we have two coupled 
linear differential second-order equations to solve: 

-h+-h--A--A--q 2 A = , (51) 

2£ 2 At 2V ' y J 

& + ifc = 2A + Jj(A-A) • (52) 
Subtracting these two equations we have: 

h = -3£A- ^A-(g 2 - ^)a + F(A) , (53) 

where 

F(X) = 4£A - ^A . (54) 



With the equation ( j53|) and equations (pT[) , we obtain the differential 
third-order non-homogeneous equation for A(£): 



where 



^(A) = ^F(A) + ijF(A) • (56) 



Performing the transformations: 

. k . 
A = — and k = £7 , 
£ 

we obtain the following differential second-order non-homogeneous equation: 

^ + k 1 + k/ 1 = Hx) ' (57) 

whose solution for A(£) is: 
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+\c 3 , (58) 

where z = 2gt^/y / 3 and W(^, — |) = J'iJ-± — J\J'-± the Wronskian con- 
structed with the Bessel functions. 

Equations ( (5"U| ) and (^) are the solutions for scalar field perturbation and 
for density perturbation. 



AW = T 



5 
£2 



3.4 Incoherent Matter (a = 0) 

Here, we have the following set of coupled differential equations: 

S + ^ = >-A)(^) + 2A + 4X^0(l +U , , ,58, 

+ ^A = ^(3 + 2w)A , (60) 
A-(ifc-fl7jf)=0 , (61) 

<5l> j + —8U j = . (62) 

Setting the four-velocity perturbation 5U k null , which, in this case, is 
allowed by one infinitesimal gauge transformation, we have: 



2r • • ilr — 2) 
(A-A) + T (A-A) + ^ 1 ^(A-A) = , (63) 
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whose solution is 

A = A + t 2(1_r) . 
Substituing (|64D in (|60|) , we obtain: 



A + -A + — q 2 X 



6r 2 - 13r + 6 



t t 2r * t 2r 
Now, perfoming the transformations: 

x = t p with p = 1 — r , A = x m g with m = - 



2(1 -r) 

we obtain the differential second-order non- homogeneous equation: 

1 



y 2 g" + yg' 



y 



2(1 -r) 



2 _ 6r 2 - 13r + 6 
3 ~ (1 _ r )2 



Its solution is given by: 

6r 2 - 13r + 6 



0(2/) 



(1-r) 2 

i / My) 



y 2 \W(u, -v 



My) 

dy 



J-u(y) 



dy+ 



+ J -Ay) I - 

J y 2 
where v = —to. 

The solutions for \(t) and A(t) are: 

6r 2 - 13r + 6 



y 2 \W(i>, -z/), 
+C 1 J„(y) + C 2 J-„(y) 



X(t) 



qt? 

+J-u(y) 



My) 



i / J-u(y) 



t 2 - r \W{u,-u) 

■My) 



dt+ 



dt 



A(f) 



t 2 - r \W(is,-v), 

+ (<V«,(i/) + <V-,,(y))H , 

6r 2 - 13r + 6 \ , , , f 1 / J_„(y) 



My) 



t 2 ~ r \W(v,-v) 



\dt+ 
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£ 2 ~ r Viy(z/,-z/) '" 
\-(t-t)[C 1 J v (y) + C 2 J-„(y)) +t 2(1 " r) • (69) 



The equations above are the solutions for A(t) and A(t) for r ^ 1. 
On the other hand for r = 1 equation (|63|)-(l65|) admits the particular 
solution A = A + C^ 1 + Cte, where 



A = r, m= -i±v^t? (70) 



4 Asymptotic Behaviour 

In this section we analyse the solutions obtained above by calculating their 
asymptotic expressions ||15||. We can readly see that the arguments of the 



Bessel functions found in the last section are indeed the ratio beetwen the 
physical wavelength Aj = a/q of the perturbation and the particle-horizon 
distance if -1 . We have, then, Xf > H^ 1 when t — > and Xf < H^ 1 when 
t — > 00. These relations invert when r > 1. 

In what follows, for sake of simplicity, we have eliminated all gauge modes 
appearing in the solutions [T^], which are represented by the constant C 3 , and 
we have chosen null phases for the oscillating terms. 

4.1 Vacuum 

i) t -»• 

A(t) ~ C^ 1 -^ + C 2 , (71) 

h{ t) ~ _^f 2 t 1 - 2r + g 2 t 5 " 4r - 2rt- 2 'V 
1 — 3r V / 
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+ T ^{q 2 t 3 - 2r + 2t- 1 + (A + 2r)r^ ; (72) 

ii) t — > oo 

A(t) ~ t 1 ^ [Ai cos^ 1 - 7 ') + A 2 sin(t 1 - r )] , (73) 
h(t) ~ t^ 1 [Ax sin(t 1 - r ) + A 2 cos(t 1 - r )] , (74) 
with Ai and A 2 given by some functions of C\ and C 2 . We also recall that 



c+l Tv /l + ^ 
4 + 3c; 



4.2 Inflation 



i) for t ^ 



A(f) -d^-^ + ^r^) , 



cos(t) 



sin(t) 



(75) 
(76) 



i=i l j i=i 

ii) for t — > oo 

A(t) ~ t- 2(1+r) [d cos^) + C 2 sin(t^ r )] , (77) 
MO-EA^ + EA^ , (78) 

i=l i=l 

where A^B^ Cj and Dj are constants and the exponents are given by: 



xi = 2 — 3r , 
X5 = 5 |r , 
yi = -1 - 4r 

1/5 
Z\ 

z 4 = -\-2r , 



l-fr 
-|-2r 



X 2 

Xq 
V2 



5r , 
7 



_1_ 13 r 

_5 _ H r 
2 2 ' 



-|-3r , 
2r , 



2 

-2 



^3 

Z3 

z 6 



|-3r , 

11 _ 11 

2 2 ' 

-3-|r , 

_I _ 15 r • 

2 2 : 



-§-2r 



X4 — 3 2 ^ ? 



1/4 



-3-fr 



= -4r 



17 



The constants above are all combinations of C\ and CV Finally we recall 
that 

1 

r = u + - . 

4.3 Radiation 

i) t 

A(t) ~ Ci + C 2 ti , (79) 
A(t) ~ C x t + C 2 ^ ; (80) 

ii) £ — >■ oo 



A(*) ~r2[Cisin(^) + C 2 cos(t5)] , (81) 

A fc (t) ~ C[ sin(t5) + C 2 cos(^) . (82) 

The non-homogeneous part of the equation ( |55|) gives, by a direct inspec- 
tion, a mode that in the asymptotic regime decays with time. 

4.4 Incoherent Matter 

i) t ^0 

A(t) ~ Ci + Cat -1 + A(r, g)£ r ~ 2 , (83) 
A(t) ~ A(t)+i 2(1 ~ r) ; (84) 

ii) t — > oo 

A(t) ~ ^ [d sin(t 1 - r ) + C 2 cos(t 1 - r )] + sin^ 1 ^') , (85) 

A(t) ~ A(t)+t 2(1 ~ r) , (86) 
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we recall that 



2 + 2u 

r = 

4 + 3a; 



5 Analyses of the Results 
5.1 Vacuum Case 

In the vacuum case we verify that the exponent of the background solutions 
is always less than unity so that, for t — > 0, both X(t) and h(t) have growing 
modes. On the other hand for t — > oo, \{t) and h{t) have only oscillating 
decaying modes. It means that perturbation of the scalar field and the per- 
turbation of gravitational field whose wavelengths are greater than particle- 
horizon are growing perturbation in contrast with those having wavelengths 
smaller than the particle-horizon which are all decaying perturbations, as it 
is an expected result. This result generalizes those obtained in [O . 



5.2 Inflation Case 

We have in the present phase three cases: r < 1, r = 1 and r > 1. (The 
case r = 1 represents the transition between two regime in the same phase: 
inflationary expansion and "soft expansion".) 

a) r = 1: this can represents the transition between two regime 
along the same phase: inflationary expansion and "soft expansion". From 
equation ([44]) we have, 

i) for q < 2, \(t) has decayind modes, 

ii) for q > 2, \(t) has oscillating modes. 
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b) r < 1: there is no inflationary inflation, in spite of the fact that 
a = — 1. Moreover, when u < —1/2, we will find an important growth of A(t), 
but a decaying scale factor. For t — > 0, \(t) has one growing mode and one 
decaying mode, while h(t) has only decaying oscillating mode. For t — > oo, 
X(t) has only decaying oscillating mode and h(t) show growing modes. 

c) r > 1: this is the case where there is strictly inflation. The 
scalar filed grows at power t 2 and the scale factor undergoes huge expansion 
for uj » 1/2. For t — > all modes of A(t) are oscillating decaying and h{t) 
has gentle grow modes for 1 < r < 2 (1/2 < uj < 3/2). For t — > oo, \(t) and 
h(t) are completly decreasing perturbations. 

5.3 Radiation Case 

In this case, r = 1/2 and the background evolves as if the scalar field were 
constant. For t — *■ 0, X(t) and A(t) (density perturbation) undergoes gentle 
growing, as t and t^. For t — > oo both perturbations show only decay- 
ing oscillating modes. In this radiation case the scalar perturbation gives a 
decaying contribuition to the corresponding GR results. 

5.4 Incoherent Matter Case 

Here, we have also three case r = 1, r < 1, andr > 1. The behaviour of A(t) 
is the same of X(t) plus one term in form of power- law. 

a) r = 1: from particular solution from equation ( |65D we see that 
i) for q < 1/2, \(t) has decaying modes; 
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ii) for q > 1/2, A(£) has two oscillating modes with decaying 

amplitude. 

A(£) has a supplementary constant term. 

b) r > 1: in this case for t — > 0, A(£) has oscillating modes, one 
of then with growing amplitude. The behaviour of A(£) shows also a pure 
decreasing mode. For t — > oo, A(£) and A(£) have only decaying modes, 
except for the case r > 2 (—3/2 < u; — 4/3). We note that in this case we 
have also a superluminal expansion of the background. 

c) r < 1: only A(£) has a growing mode for £ — > 0. A(£) has only 
decaying oscillating modes and A(£) is still increasing, for £ — > oo. In the case 
of —4/3 < to < — 1 we have r < and then A(£) is growing with high power 
of £. We observe that, for this case, we have a contracting Universe. So it 
can be consider of physical interest if it can be applied to some primordial 
phase. 

Our results confirm those obtain in fl2~0|), where it is pointed out that the 
BD solution of this case can not leads to a solution of the structure formation 
problem, when we consider values of u imposed by local physics. 

6 Conclusion 

We gave in this paper a general classification of perturbation in BD Theory 
corresponding to the phases of the expanding universe where the matter is 
described by equations of state p = ap with a = 0, a = —1, a = 1/3 
and p = 0. We assumed that the background solution for the scalar field 
and for the scale factor are that of "pure" power-law. We have found the 
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exact solutions of the perturbed equations in the form of integrals of Bessel 
functions. The calculation here confirm, in general, the very known fact: 
perturbations at scales larger than particle horizon distance grow up, while 
for smaller scales the perturbation exhibits an oscillating behaviour. It is 
not true for incoherent matter where density perturbations can be always 
amplificated. 

In the inflationary phase we can see that the strong expansion of the 
universe erases the perturbation on the scalar field but the gravitational 
perturbation is growing for scales greater than particle-horizon distance. On 
the other hand, it is very important to emphasize that the perturbations in 
the density of the matter are null in Brans-Dicke inflation. This recovers the 
main results from classical inflationary scenario. 

The radiation phase do not give an expressive contribution with respect to 
the GR corresponding results. In the incoherent matter the situation is more 
rich. For some negative values of uo we can have a superluminal expansion 
and at the same time a significant amplification of perturbations. 
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